Density Functional Theory Characterization of the Multiferroicity 
in Spin Spiral Chain Cuprates 
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The ferroelectricity of the spiral magnets LiCu2C>2 and LiCuVCU was examined by calculating 
the electric polarizations of their spin spiral states on the basis of density functional theory with 
spin-orbit coupling. Our work unambiguously reveals that spin-orbit coupling is responsible for the 
ferroelectricity with the primary contribution from the spin-orbit coupling on the Cu sites, but the 
asymmetric density distribution responsible for the electric polarization occurs mainly around the 
atoms. The electric polarization is calculated to be much greater for the ab- than for the bc-plane 
spin spiral. The observed spin-spiral plane is found to be consistent with the observed direction of 
the electric polarization for LiCuVCU, but inconsistent for LiCu2 02- 

PACS numbers: 75.80.+q,77.80.-e,71.15.Rf,71.20.-b 



In the past several years there has been a revival of in- 
terest in understanding magnetic ferroelectrics because of 
their potential applications in novel magnetoelectric and 
magneto-optical devices [H, H, B 0, II] • In these multifer- 
roics, the paramagnetic phase is centrosymmetric, and 
electrical polarization appears only at the transition to a 
magnetically ordered phase, which is responsible for re- 
moving the inversion symmetry thereby generating a po- 
lar field. Two mechanisms of multiferroicity have been 
proposed in the literature. The exchange striction due to 
the symmetric parts of the exchange coupling appears rel- 
evant for the multiferroicity in systems such as RM^Os 
Q . The ferroelectricity (FE) in spiral magnets has been 
suggested to originate from spin-orbit coupling (SOC) 
0,1a H EH- However, this suggestion has never been 
verified by first principles electronic structure studies. 

Very recently, FE was discovered in two cuprates 
LiCu2C>2 [ll| and LiCuVC>4 These oxides contain 
spin-frustrated CuC>2 ribbon chains that are made up 
of edge-sharing CuC>4 squares with spin-i magnetic ions 
Cu 2+ , for which an unpaired spin resides in the e g or- 
bital (commonly referred to as the d x 2_ y 2 orbital, but 
the d xy orbital in the local axis system adopted in the 
present work). Due to the competition between nearest- 
neighbor (NN) ferromagnetic and next-nearest-neighbor 
antiferromagnetic interactions in each CuC>2 chain, a spin 
spiral state can set in along the chain direction (i.e., the 
b direction) at low temperature [l3| . Neutron diffraction 
studies showed that the Cu 2+ moments lie in the CuC>2 
ribbon plane (i.e., the afe-plane) in th e sp in spiral state 
of both LiCu 2 2 H and LiCuV0 4 15]. Results of a 
recent ESR study [16| on LiCu 2 2 are consistent with 
the a6-plane spin spiral. Below the temperature of the 
spiral-magnetic order, the electric polarization is found 
to be along the a direction in LiCuVC>4 [12j, but along 
the c direction in LiCu2C>2 [Hj]. This difference in the 
electric polarization directions is puzzling because, ac- 
cording to the Katsura-Nagaosa-Balatsky (KNB) model 
, the electric polarization along the a-direction is asso- 



TABLE I: Electric polarizations (in units of (iC/m 2 ) of 
LiCu202 and LiCuVCU calculated for the ab-, ac-, and bc- 
plane spin spiral states of Fig.[2]using the experimental (EXP) 
centrosymmetric structures and the optimized (OPT) struc- 
tures. 

P(ab) P(ac) P(frc) 
LiCu 2 2 (EXP) (51.4,0,0) (0,0,0) (0,0,-9.8) 
LiCu 2 2 (OPT) (441.4,0,0) (0,0,-191.8) 
LiCuV0 4 (EXP) (103.5,0,0) (0,0,0) (0,0,-15.7) 
LiCu2Q4 (OPT) (595.6,0,0) (0,0,-223.0) 
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-97.18 




FIG. 1: Total energy E(q) calculated for the spin spiral state 
of LiCu202 as a function of q = (0.5,^,0) on the basis of 
the non-collinear LDA+U method using a unit cell containing 
four formula units. 



ciated with the etfr-plane spin spiral, and that along the 
c-direction with the fee-plane spin spiral. However, it is 
unclear whether the KNB model is applicable to the e g 
systems [l2} because it was formulated for a t 2g system. 
Furthermore, in their recent s tudy of the spiral magnet 
LiCuV04, Jia et al. reported 17[ that the SOC on the 
ligand O 2p orbitals plays a more important role than 
does the SOC on the Cu 3d orbitals. It is important to 
verify whether or not this conclusion is correct. 
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FIG. 2: (color online) Spin arrangements in the (a) ab-, (b) 
ac-, and (c) 6c-plane spin spiral states using the 1x4x1 
supercell. The large white and the small red circles represent 
the Cu and O atoms, respectively. The legends and ® 
denote z and — z directions, respectively. The directions of 
the electric polarization corresponding to a given spiral plane 
are also shown. In the present work, the a, b and c axes 
correspond to the local x, y and z axes, respectively. 



In this Letter, we address the aforementioned issues on 
the basis of non-collinear fully relativistic density func- 
tional calculations for LiCu2C>2 and LiCuVC>4. The elec- 
tric polarizations associated with the non-collinear spin- 
orbit coupled states are calculated using the modern po- 
larization theory [l8| for the first time. To the best of 
our knowledge, this is the first study on multiferroicity 
in spiral magnets using first principles method with the 
SOC effect included. Our electronic structure calcula- 
tions were carried out using the local density approxima- 
tion (LDA) on the basis of the projector augmented wave 
method [l9| encoded in the Vienna ab initio simulation 
package (VASP) [13] ■ The plane- wave cutoff energy was 
set to 400 eV. To properly describe the strong electron 
correlation in the 3d transition- metal oxide, the LDA plus 



on-site repulsion U method (LDA+U) was employed [21 1 . 
In the following, we report results obtained with U = 6 
eV and J = eV on Cu, but our results with other U 
values between 4 — 7 eV are similar. 

In the absence of SOC, the incommensurate spiral 
magnetic order can be simulated without resorting to 
the supercell technique due to the generalized Bloch the- 
orem 22]. A neutron diffraction study [3] showed that 
the spin structure of LiCu202 is helimagnetic with prop- 
agation vector q = (0.5, 0.174, 0). We performed a series 
of calculations [23J for the spin spiral state with prop- 
agation vector q = (0.5,^,0). The dependence of the 
electronic energy E(q) upon q y is presented in Fig. [TJ 
which shows that the energy at q y — (i.e., the FM 
arrangement between all NN's) is lower in energy than 
that at q y ~ 0.5 (i.e., the AFM arrangement between 
all NN's). This is in accord with the FM nature of the 
NN exchange interaction. The energy minimum occurs 
at around q y = 0.165, in good agreement with the exper- 
imental value of q y = 0.174. This agreement indicates 



that the spin spiral ground state of a system with spin 
frustrated Cu02 ribbon chains can be predicted by non- 
collinear LDA+U calculations. 

To see if the LDA+U method can describe the FE in 
LiCu202 as well, we calculated the electric polarization 
of the spin spiral state with q y = 0.174 using the Berry 
phase method [la ]. This calculation leads to negligible 
electric polarization. As a possible reason for this fail- 
ure, we considered the lack of geometry relaxation in 
the spin spiral ground state. However, LDA+U calcu- 
lations for the fully optimized structure of LiCu202 does 
not lead to any appreciable electric polarization. Simi- 
larly, our LDA+U calculations for the spin spiral state 
of LiCuV04 show no electric polarization. Consequently, 
the exchange striction is not an appropriate mechanism 
of electric polarization for LiCu202 and LiCuVOit. 

The above observation led us to examine SOC effects 
on the electric polarization in the spin spiral states of 
LiCu 2 2 and LiCuV0 4 . We carry out LDA+U+SOC 



calculations for the q = (0,0.25,0) spiral states [24 1 
with three different spin spiral arrangements, shown in 
Fig. [2j by using a 1 x 4 x 1 supercell. Our calculations 
for LiCu2 02 give rise to substantial electric polarizations 
(see Table [U). The electric polarization P is calculated to 
be along the a-direction for the afr-plane spin spiral, and 
along the — c-direction for the frc-plane spin spiral [31 1. 
When the spin spiral is in the ac-plane, the electric po- 
larization is calculated to be zero. The same results are 
also obtained for LiCuV04. Thus, for both LiCu202 and 
LiCuV04, the directions of the calculated electric polar- 
izations are consistent with the prediction of the KNB 
model. This agreement is somewhat surprising since the 
KNB model was derived for a linear trimer M-O-M with 
t2 g transition metal ions M jjjj]. The calculated electric 
polarizations are anisotropic, e.g., the electric polariza- 
tions of LiCu202 are 51.4 and 9.8 [iCjrr? for the ab- 
and 6c-plane spin spiral arrangements, respectively. Note 
that the anisotropy in electric polarizations is not pre- 
dicted by the KNB model. 

Our calculations of electric polarizations show that the 
observed spin-spiral plane is consistent with the observed 
direction of the electric polarization for LiCuV04, but 
this is not the case for LiCu2 02- It is of interest to 
consider a probable reason for the latter discrepancy. 
For LiCuV04 (2||, the planes of their spin spiral can be 
flipped by the application of an external magnetic field. 
For LiCu202 the direction of the electric polariza- 
tion can be flipped by an external magnetic field. In 
the ideal crystal structure of LiCu202, the Cu02 ribbon 
chains with Cu 2+ ions are interconnected by linear O- 
Cu-0 bridges with diamagnetic Cu + ions. Matsuda et al. 
reported that the actual composition of their "LiCu202" 
sample is given by Li1.i6Cu1.s4O2.01, in which 16% Cu 2+ 
sites of the Cu02 ribbon chains are replaced with dia- 
magnetic Li + ions due to a good match of the ionic radii 
of Li + and Cu 2+ while 16% Cu + ions in the linear O-Cu- 
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FIG. 3: (color online) (a) Cross section view of the electron 
density difference between the LDA+U+SOC and LDA+U 
calculations for the afe-plane spin spiral state (Fig. HJa)) 
of LiCuVCU. (b) Perspective view of an isosurface cal- 
culated for the electron density difference between the 
LDA+U+SOC and LDA+U results for the bc-plane spin spi- 
ral state (Fig. [2^c)) of LiCuV04. The red and blue surfaces 
represent —0.25 and 0.25, respectively. 

O bridges become Cu 2+ ions due to the charge balance 
requirement. This implies that the spins of the Cu 2+ ions 
residing in between the CuC>2 ribbon chains can act as a 
source of "external" magnetic field for the CuC>2 ribbon 
chains, and hence can influence the plane of spin spiral 
in the CuC>2 ribbon chains. Similarly, the Ising behav- 
ior of CuFe02 with high spin Fe 3+ (d 5 ) ions below 14 K 
has been explained in terms of a uniaxial magnetic field 
associated with oxygen defects [261 ] - Since the nonstoi- 
chiometry of "LiC^CV samples would depend strongly 
on synthetic conditions, it is highly desirable to measure 
the plane of spin spiral and the direction of electric po- 
larization by using the same sample. 

To test the conclusion of Jia et al. [17| that the SOC on 
the O 2p orbitals is more important for the ferroelectric 
polarization of LiCuVC^ than is the SOC on the Cu 3d 
orbitals, we carried out LDA+U+SOC calculations for 
LiCuV04 with the afr-plane spin spiral state by switch- 
ing off the SOC on either the Cu or the O sites. When 
the SOC is neglected on the Cu sites but kept on the O 
sites, the electric polarization is calculated to be along 
the -a axis with P a = -29.8 [iC/m 2 . If the SOC is kept 
on the Cu sites but neglected on the O sites, the elec- 
tric polarization is calculated to be along the a axis with 
P a = 127.1 fiC/m 2 . The sum of the above two polariza- 
tions is very close to the electric polarization (P a — 103.5 
/iC/m 2 ) calculated when the SOC is kept on both Cu and 
O sites. Therefore, the primary contribution to the elec- 
tric polarization comes from the SOC on the Cu sites, 
which is greater than that from the SOC on the O sites 
by a factor of approximately four. 

The dominance of the SOC on the Cu sites does not 
necessarily mean that the asymmetric charge distribu- 
tion, needed for nonzero electric polarization, is associ- 
ated mainly with the Cu atoms. To examine how the 
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FIG. 4: (color online) Difference in the PDOS plots of the 
adjacent Cu atoms, i.e., the Cu2 and Cul atoms of Fig. [2ja), 
in the (a) ab- and (b) &c-plane spin spiral state of LiCuV04. 
The energy region plotted corresponds to the hole xy state of 
LiCuVCU. 



FE arises from the spiral magnetic state, we analyze the 
electron density distribution of the spin spiral state of 
LiCuV04 by plotting the difference between the electron 
density of the LDA+U+SOC calculation and that of the 
LDA+U calculation. For the case of the afr-plane spin 
spiral, Fig. 02a) shows that the difference density around 
each Cu ion is almost symmetric, and hence contributes 
little to the electric polarization. However, an asymmet- 
ric difference density is found on the O atoms such that a 
nonzero electric polarization sets in along the a-direction. 
Fig. &h) shows a perspective view of the difference den- 
sity plot for the case of the 6c-plane spin spiral. Again, 
the difference density around each Cu atom is nearly sym- 
metric, and an asymmetric difference density is found at 
each O atom largely along the c-direction eventually lead- 
ing to the electric polarization along the c-direction. 

From the viewpoint of local bonding, the SOC at each 
Cu site mixes into the hole xy state other 3d-states (i.e., 
the x 2 — y 2 7 xz, and yz states). The extent of the mix- 
ing for each 3d-state depends on the spin direction at a 
given Cu site because of the SOC term XS ■ L. This gives 
rise eventually to the asymmetric electron density distri- 
bution on the O sites because each 3d state of a given 
Cu site has the 2p orbitals of its surrounding O atoms 
combined out-of-phase with the Cu 3d orbital. We con- 
firmed this point in two ways. First, we analyzed the 
partial density of states (PDOS) for the hole states of 
LiCuV04 (i.e., the unoccupied band associated with the 
down-spin xy state from each Cu site), which are located 
in the energy range from 1.1 to 1.8 eV above the Fermi 
level. In the spin spiral state, the neighboring Cu ions in 
each Cu02 chain with different spin directions will have 
slightly different 3d PDOS because the mixing between 
the 3d-states induced by the SOC depends on the spin 
direction. As expected, the Cul and Cu2 ions [as defined 
in Fig. [2Ja)] show differences in their PDOS not only for 
the ab- plane spin spiral (Fig. EJa)) but also for the fee- 
plane spin spiral (Fig. HJb)). Second, we extended the 
tight-binding calculations of Jia et al. Q to a tetramer 
model in which two Cu ions are bridged by two O ions 
as found in the Cu02 ribbon chains, and carried out nu- 



4 



merical solutions by using all 3d orbitals of Cu and all 
2p orbitals of O. As expected, the occupations of the 
x 2 — y 2 , xz, and yz states at a given Cu site are found 
to depend on the spin direction of that site. In addition 
to a transverse electric dipole moment, a nonzero polar- 
ization is found along the Cu chain direction. This latter 
longitudinal component oscillates as the spin spiral pro- 
gresses along the chain, and hence its effect is canceled 
out. Consequently, the overall polarization is calculated 
to be along the a and c directions in the case of the ab- 
and fee-plane spin spiral arrangements, respectively. 

Concerning the effect of SOC on the FE in spiral mag- 
nets, there is a debate as to whether the FE arises from 
pure electronic effects @, [E_l or displacements of ions [_| . 
Using first principles calculations without taking SOC 
into account, Picozzi el al. [___| showed that both factors 
are active for the FE in HoMnOs. In our calculations de- 
scribed above, we used the ex per imental centrosymmetric 
structure for both LiCu 2 2 28] and LiCuV0 4 [___|. To 
evaluate the effect of polar atomic displacements on the 
electric polarization, we optimized the crystal structures 
by LDA+U+SOC calculations until the atomic forces be- 
come less than 0.01 eV/A[32j]. The calculated electric 
polarizations using the resulting relaxed structures are 
shown in Table HJ which shows that the direction of the 
electric polarizations remains unchanged but the mag- 
nitude of the electric polarization becomes strongly en- 
hanced by the structural relaxations. It is found that the 
Cu 2+ ions move along the —a and c directions in the case 
of the ab- (Fig. _2a)) and fee-plane (Fig.^c)) spin spiral 
states, respectively. With respect to the centrosymmetric 
structure, however, the largest displacement of the Cu 2+ 
ions is about 0.0004 A, which are hardly detectable in 
experiment. Note that the experimental electric polar- 
ization at 2.5 K for LiCu202 is 4 /j,C/m 2 , which is rea- 
sonably well reproduced by the calculation using the ex- 
perimental structure (i.e., 9.8 fiC/m 2 ). The electric po- 
larization 191.8 fiC/m 2 calculated by using the relaxed 
structure is much greater. Since the polar atomic dis- 
placements associated with the relaxation are extremely 
small, it is probable that their effects are negated by zero 
point vibrational effects. 

To probe how one might improve the effect of multi- 
ferroicity, it is necessary to examine how SOC affects 
electric polarization. For this purpose, we performed 
LDA+U+SOC calculations for LiCuV04 by treating the 
speed of light c as an adjustable parameter and noting 
that the strength of SOC is inversely proportional to the 
square of the speed of light c These hypothetical 

calculations show that the magnitude of electric polar- 
ization is proportional to the strength of SOC. Since the 
SOC of an atom increases in strength with increasing 
the atomic number, high-temperature multiferroics with 
large electric polarization are expected from spin-spiral 
systems with 4d or 5d transition metal elements. An- 
other important factor is that the critical temperature of 



the spin spiral ordering should be high. 

In summary, our study of LiCuV04 and LiCu202 re- 
veals that SOC is responsible for their FE, the primary 
contribution to the FE arises from the SOC on the Cu 
sites, but the asymmetric density distribution responsi- 
ble for the electric polarization occurs mainly around the 
O atoms. For both LiCuV04 and LiCu202, the electric 
polarization is calculated to be much greater for the ab- 
than for the fee-plane spin spiral. The observed spin-spiral 
plane is consistent with the observed direction of the elec- 
tric polarization for LiCuV04. However, this is not the 
case for LiCu202- To resolve the latter discrepancy, fur- 
ther experimental and theoretical studies are necessary. 
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